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I. INTRODUCTION 


The subcritical assembly initially was utilized as an 
experimental unit to determine if a chain-reacting assembly 
were possible. As reactor technology progressed the assembly 
came to be used more and more as a test unit for proposed 
reactor designs and materials, and the majority of these test 
investigations have been conducted in the steady or equili- 
brium state. Operation of the assembly by disturbance of the 
steady state has led to new areas of investigation, 

Frequent utilization of the subcritical assembly as a test 
unit has come about naturally. Cost of a subcritical assembly 
is considerably lower than that for a critical assembly since 
special safety equipment is not needed (11, p. 103). Less 
fuel and moderator are used, and experiments can be conducted 
easier in the subcritical than in the critical assembly. 

Disturbance of the steady state may be accomplished by 
two methods, reactivity forcing and source forcing. Re- 
activity forcing can be represented as a change in the neutron 
Gensity brought about by a small harmonic variation in the 
reactivity, whereas source forcing is produced by a variation 
in the rate at which neutrons are being added to the assembly. 

Thus far the investigation of the dynamic response of a 
system has been limited to servomechanism theory, which 


treats the system as a “black box". Phase and amplitude 
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response of the system in relation to a forcing function must 
be determined in order to apply servomechanism theory. The 
above relationship called the transfer function is derived by 
noting the responses as indicated brought about by either 
reactivity or source forcing. 

it is believed however that neutron diffusion in a sub- 
critical assembly will not produce identical phase and 
amplitude responses throughout. These variations should depend 
upon the method of varying the rate of neutrons from the 
source and the position of neutron-flux measurement. 

The purpose of this investigation was threefold. Firat, 
an experimental verification of tne theoretical transfer 
function for a sudcritical assembly was attempted using a 
varying rate of neutrons from the source. Secondly, neutron 
diffusion in the subcritical assembly that shows variance 
from the predicted identical phase and amplitude response 
throughout the assembly was investigated. The last area of 
investigation was an attempt to determine assembly character- 
istics by using a sinusoidally varying source. 
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Il. REVI8W OF THE LITERATURE 


A reasonable amount of literature is available on the 
development and usage of the transfer function of a critical 
assembly. Original work in the field was carried out by 
Qwens, Crever and Pigott (12) in 1945, in which an automatic 
control unit was used to compensate for reactivity variations 
in order to maintain the desired power ievel in a reactor. 
Their work showed that the assembly transfer function ex- 
hibited proper characteristics to deseribe the reactivity 
variations in the assembly. 

Later that same year, Franz (4) developed the general 
form of the trenefer function of a nuclear reactor by apply- 
ing the standard techniques of servo theory. He applied the 
transfer function to an electronic assembly simulator te pre- 
dict assembly response to a step-function chance of 
reactivity. 

The practical verification of using a subcritical assembly 
in testing reactor components was shown by Axtmann, Dessauer, 
and Parkinson (1) in 1955. A chain-reacting "test assembly" 
was utilised to tert materials by noting effects on the Kinet- 
ie behavior of the assembly. Identical tests were conducted 
in a subcritical assembly, and the results were nearly 
identical with the subcritical assembly producing results 
faster, cheaper and safer. 
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Glasstone (5) briefly discussed the transfer function in 
a subcritical assembly. In his development he assumed a 
sinusoidal variation in source strength and used electrical 
analogs of reactor quantities to indicate the relationship of 
electrical parameters, with the reactor paraueters, reactivity 
forcing and source forcing. 

Weinberg and Schweinler (15) developed the equations for 
the response of a critical chain~-reacting assembly to a 
thermal neutron absorber which is oscillated back and forth 
ineide the assembly. They further stated that at frequencies 
which are low compared to the reciprocal periods of the de- 
layed neutrons, the neutron density in the assembly rises 
ami falls as a whole. However, if the frequency is much 
greater than the reciprocal period of the delayed neutrons, 

@& spherical neutron wave emanates from the vicinity of the 
oscillator and is attenuated in an exponential manner with 
distance as 1t is propagated. 

Weinberg and Wigner (16) extended this theory to sub- 
critical assemblies. It was stated that the neutron-wave 
wave length is related to the material buckling. They sug- 
gested that the exponential~experiment method of determining 
material buckling could be performed by measuring the proper~ 
ties of neutron waves established by a localized oscillating 
neutron source. A particularly desirable feature 
of this method is that an experiment to determine material 
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buckling can be performed in an assembly of smaller dimensions 
than that required by the normal exponential experiment. 
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Tit. LIST OF SYMBOLS 


Activity 

Width of assembly 
Attenuation length 
Buckling 


Fraction of average group of delayed 
neutrons 


Fraction or 15> group of delayed neutrons 


Concentration of average group of delayec 


neutron precursors 
End-correcticn factor 
Hermonic«correction factor 


Concentration of 1°" group of delayed 
neutron precursors 


Thermal neutron diffusion coefficient 
Thermal utilization 
inverse relaxation length 
\=a 
Constant 
Effective multiplication factor 


Multiplication factor for an infinite 
assembly 


Neutron diffusion length in assembly 


Average neutron lifetine in a Pinite 
assembly 


Mean neutron lifetime in an infinite 
assembly 


countsa/min, 
in. 
in., cm. 


em. 7* 


Cit. 


th. ~* 


Cm. 


Bec. 





A Decay constant of average group of de- 
layed neutron precursors 
Aq Decay constant of the 4 bh group of de- 
layed neutron precursors 
M= Migration area 
n Number of neutrons of thermal energies 
(J Frequency of oscillation 
p Resonance escape probability 
g Thermal neutron flux 
a Macroscopic cross section 
t Time 
Ta Fermi age 
Vv Thermal neutron mean velocity 
Zz Vertical distance above base of assembly 
Subseript: 
& Absorption 
Gn Delayed neutron 
e External 
& Graphite 
i qth part 
m Harmonic mn 
p Prompt 
Q oteady state 
i Complex amplitude 


11,13531,33 


Harmonics of Y 


sec.7* 


sec.7* 
in.*, om,* 
neutrons /cm.? 


radians 
sec. 


neutrons 
cm." sec. 


in,**, em, ** 
BeC. 
om. * 
em./sec. 


in. 
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IV. THEORETICAL ANALYSIS 
A. Transfer Function 


The transfer function is defined as the ratio of the 
phase and amplitude response to a sinusolidaliy-varying forcing 
function. This response for a subcritical assembly can be 
obtained from the time-dependent general diffusion equation. 
Thermal neutrons only will be considered since the experi- 
mental equipment is designed to vary the rate at which thermal 
neutrons enter the assembly. The basic equation is (6, 


yp, 101) 





‘ a | on 1 2 2 7 
where 7 = nv, and v is assumed to be constant. 

The source term can be separated into three parts. The 
first term is due to the prompt neutrons resulting from 
fission. Fermi-age theory yields this part of the source 
term to be (5, p. 226) 

=p* 7 
where G is the fraction of all the delayed neutrons. 

The second part of the source term is due to the delayed 
neutrons resulting from decay of various fission products. 
A total of six groups of delayed fission neutrons are 


utilized in expressing this part of the souree term. Each 
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group has its own value of hes 6B, end C4, The net rate of 


formation of the delayed-neutron precursers is then 


K 60 








O° al 6, 


— Qo B= Xz Cy . : Eq. 3 


The delayed-neutron source term using Fermi-age theory can be 
expressed as (5, p. 227) 
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where e Aq C4 is the decay rate of the precursor which is 
Lm 
numerically equal to the rate of production of the delayed 
neutrons. The Fermi-age for the delayed neutrons has been 
taken equal to that of the prompt neutrona (5, p. 227). 
The last part of the source tezm is due to the external 
source which in this analysis supplles the source forcing 
needed. it is written as 
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With substitution of the sbove expressions of the source 


term into Equation 1 the following equation is obtained 
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The space and time variables are considered separable 
when kK remains constant (5, p. 227). This yields an ordinary 
differential equation for Equation 6. If the subcritical 
system is large, then v* @ = -B® @ (6, vo. 361). 
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and divide Equation 6 through by 2,, to obtain the result 
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With division of Equation 9 through by (1 + L®B*) and the use 
ef the relationships 
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With division by /, and utilization of @ =nv, Equation 
i2 is reduced to 
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With rearrangement Equations 1¢ and 11 became 
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In a sinilar manner, Equation 3 with the use of Equations 


10 and 15 is found to be 
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Solving for ) 4 C, and substituting into Equation 16 yields 
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In this particular case the value of k will remain 
constant. Neither assembly geometry nor fuel-loading changes 
will take place in the experiment. However both the neutron 


density and precursor concentration will vary as 
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nen, +n, eWV* Eq. 21 


and 


Cy =o, +0, 9 Y* Eq. 22 


where n, and C, are the complex amplitudes. With the aid of 
Equations 17, 21 and 22, Equation 20 yields 
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The complete solution to Equation 23 involves using the 
six groups of delayed neutrons each with its discrete value 
of Aas Gi» and Cy. The solution using an average group of 
delayed neutrons will be shown. 

Equation 23 can be separated into steady and sinusoidal 
components, and each constitutes a valid equality. The 
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transfer function uses only the sinusoidal component, which 


with group notation eliminated, is 
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Again usins Equation 17 ond substituting Equations 21 
amd 22 gives 
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Sy definition, the amplitude portion of the transfer 
function is the ratio of the response to the sinusoidally - 


varying Griving source. The ratio can be expressed as 
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Rearranging Equation 30 by putting the right side denominator 
under the common factor (A+ jw) (l+k), multiplying through 
by 2 + JW, then dividing by ) and simplifying gives the 


desired form 
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B. Attenuation Length of Neutron Wave 


When the rate at which neutrons entering a subcritical 
assembly is oscillated in a sinusoidal manner at a frequency 
considerably greater than the reciprocal period of the de- 
layed neutrons a neutron wave is propagated. With an average 
delayed neutron period of approximately 10 seconds, fre- 
quencies greater than one radian/second should produce these 
neutron waves. Frequencies considerably less than the 
delayed-neutron reciprocal period allow the delayed neutrons 
to remain in equilibrium with the prompt neutrons (14, 

p. 100). However, at higher frequencies the assembly be- 
haves as if there were no delayed neutrons and the delayed 
neutrons merely contribute to the background. Henee the 
neutron wave has characteristics that are dependent upon 

the "prompt" response of the assembly. These characteristics, 
wave length, velocity, and attenuation length are all 
dependent upon the frequency of oscillation. 

Weinberg and Wigner (16, p. 437) developed equations 


for the three neutron-wave characteristics. The equation 
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for the characteristic that is of importance here, attenua- 


tion length, is 


@a.il. 2 
Y 2m Bq. 32 
j Uk ~ 1)* + ( W £,)2 178 - (Ky - 1)? 1/2 


yvhere W is greater than one radian/secomi. Attenuation length 
is defined as the distance for wave amplitude to Gecrease by 
a factor of ¢. 

If a neutron wave were propagated by using a high- 
frequency sinusoidal source oscillation then attenuation 
length coulé be determined. Amplitudes of the sine wave 
in selected vertical positions could be measured, and a value 
for attenuation length would then be determined from the 
equation A = K e7B/Bei. | where A is the amplitude at Z, and 
K i8 a constant. A semi-~log plot of amplitude versus verti- 
cal position will yield a siope of 1/a.l. This value of 
1/a.l. is the "inverse relaxation length" of the complex 
buckling for the frequency investigated. If the phase angle 
is small at a given frequency, the attenuation length at that 
frequency would be approximately equal to relaxation length 
usec to calculate material buckling. Therefore, the complex 
buckling determined by this value of atteruation length 
would be approximately equal to the material buckling. 


~~ ee wm 


> «se 
- 4 ry 


““ te ee | 


Oo SS Ss Oe oF Ve Clee, Wt WY oe 
ree) Sie fee ©) eet ee @ ences fl 
+ mony 
et hee ee ee te eer + 4 = 
EL LP ee AAA Remy 
— ae a oe = see ore 
wt) = fee © fee Gee! Geet ey 
| mle ce! Srey mite. | 
7 .) () hie oe 4) ee tes 
iis a et) 1) 
lo ee i eB Be) 
ee teins coe ow ak 





































Ae ll at a i ed Jom i, 
aes oe — 4 : 
—_ oo —S 9 





| elms Ne , eee of 


- ==—=<——« <=> = aap 








1é 


V. EXPERIMENTAL EQUIPMENT 
A. Suberitical Assembly 


Tne subcritical assembly used for the experimental in- 
vestigation is shown in Figure 1 with the east cover removed. 
Cylindrical rods of graphite measuring 7 inches in diameter 
were machined to a square from 6 inches across, leaving 
rounded corners. These blocks of AGKR-grade graphite were 
then stacked in ten columns nine rows high to form the lower 
section of the assembly. Fox the top five rows cylindrical 
graphite rode 6-3/8 inches in diameter were machined to a 5 
inch by G inch cross section with rounded corners. The 
space available between the corners of adjoining blocks 
was utilized for insertion of fuel elements or measuring 
equipment. 

The graphite blocks were assembled upon a plywood sheet 
supported ty a wooden pedestal. As shown in Figure i, the 
pedestal was diviced into three accessible spaces approxi- 
mately one foot high. The two outermost spaces contained 
open-topped aluminum tanks, filled with water, which ex- 
tended the width of the assembly. They served to moderate 
the neutrons and to reduce radiations coming directly from 
the sources located in the oscillator unit beneath the 
assembly. Figure 2 shows the oscillator unit and selective- 


count device which was placed in the center access section 
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Fissure 2. The oaclilator unit and selective-~ 
count device 
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beneath the assembly. 

All four sides enc top of the assembled graphite were 
covered with 10-mil cadmiun samiwiched between a ©.375-ineh 
thickness of plywood end 0.125 inches of masonite. A4&pproxi-~- 
mately 95 pereent of all thermal neutrons were stopped ty 
this cadmium sheet, giving a “black boundary" to thermal 
neutrons. A wooden framework braced the north and south 
sides of the assembly and both the east and west faces were 
removeable for easy access. 

Natural uranium in the form of cylindrical rods 1.0 
inches in diameter and 6.0 inches long served as the fuel. 
These rods were encased by 2S aluminum cans with a 0.040- 
inch wall thickness and end caps 0.200 inches thick. Each 
¢an of uranium was then helically wound with a spacing wire 
and placed into a 61S aluminum tube, 62 inches long with an 
outside diameter of 1.375 inches and « wall thickness of 
9.535 inches. Approximately 10 feet of the 25-aluminun 
Spacing wire was needed for each tube. A total of seven 
canned uranium rods were placed in each tube to constitute 
one fuel element. With the fuel elements pleced in every 
ether hole as shown in Figure 1 an &.5=<inch square lattice 
in the lower region of the assembly was formed. This lattice 
configuration waa maintained throughout the experimental 
investigation. 
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B. Oseillator Unit 


The oscillator unit and selective-ctount device 1s show 
in Figure 2. The drun-type assemtly at the right ¢nd is 
the oscilintor unit itself and ig shown schematicaliiy in 
Fisure %. This unit was designed to contain the five one- 

urle PueBe sources and to vary the rate at which thermal 
neutrons were emitted from the source in 2a sinusoidal 
manner. 

Each cylindrical scurce measuring 1 incn in dlemeter 
amd 1-3/5 inches high was placed in the source retainer. 

A nonerotating-paraffin jacket surrounding the retainer was 
constructed of 1/1G-inch aluminum sheet rolled into 2 
Giameter of 7-3/5 inches. End capa of 1/é-inch alwainun 
were attached, giving an overall length of 14.5 Inches. 
Twenty-five pownis of white paraffin were poured into the 
jacket through a filler cap, and precautions were taken to 
reduce the possibility of shrinkage voids. Moderation to 
supply the thermal neutrons desired was provided by the 
three inches of persffin. The paraffin jacket ami the bear~ 
ings for the pattern cylinder were beth supported by the 
eource retainer which extenced to the support brackets. 

An outer rotating cylinder, the pattern cylinder, was 
@lso formed from 1/S-inch ealwainum sheet. Its dimensions 
were an outside diameter of & inehes and a length of 15 
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Figure 3. Oscillator unit and cadmium-pattern detail 
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inches. End plates made from 1/4-inch aluminum were fitted 
and screwed to a collar fitted on its bearing ocutererun 
casing. The immer-run casing of the hearing rested on the 
gource retainer for support. To give easy access to the 
inner section of the pattern cylinder the end plates were 
attached by screws. 

A gine pattern with the dimensions as shown in Figure 3 
was cul fron lO-mi1 cadmium. Three thicknesses, giving a 
total of 0.030 inches, were cemented to the immer surface of 
the pattern cylinder. Most of the thermal neutrons were 
stopped and a small fraction of the fast attenuated by this 
thickness of cadmiun (10, p. 12). 

Point A of the pattern was the wide~open position, allon- 
ing the maximum number of thermal neutrons to escape from the 
oscillator unit. The closed position, point B, stopped most 
of the thermal neutrons. When point A was at ton-dead 
center, the oscillator unit was at the sero-degree position. 
The 160<Jegree position occurred when point B was at top-dead 
eenter. Hence the exact position of the sine pattern of 
the oscillator unit was designated by a siven angle from oO 
to 360 oscillator degrees. 

Measurements taxen within 1.5 inches of the oscillator 
unit at the two extreme nositions deserlbed above showed 
that the peraffin jacket moderated ea satisfactory percentage 
of n@utrons to be used in the experimental investigation. 
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The value obtained at position A minus that at B was 5% per- 
cent of the value at position A. 


G. Counting Apparatus 


In order to measure the neutron flux during a designated 
number of oscillator-unit degrees a selective-count device 
was constructed. A positive, non-slip linkage had to be 
supplied from the oscillator unit to a position removed from 
the neutron sources. The equipment needed to select the 
number of oscillator degrees in which counting would take 
place was located at this "removed" position. 

A bicycle sprocket was bolted to the end plate of the 
pattern cylinder. A 110-inch bicycle chain was used to 
connect this sprocket on the cscillator unit to a sprocket 
of the same diameter mounted on a countershaft at the removed 
position. Position by degrees of the oscillator unit was 
then transmitted unchanged to the "removed" position where 
work could be performed without the intense radiation. 

The countershaft was a polishing-head spindle which was 
mounted on a 3/4-inch plywood extension connecting the 
oscillator unit to the remainder of the equipment. At one 
end of the spindle shaft a Geinch V-belt pulley was locked 
to the shaft, and the bicycle sprocket was in turn bolted 
to this pulley. At the opposite end of the shaft was locked 
@ plywood plate, 6-3/4 inches in diameter. The polishing- 
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head spindle was driven by a 1/4-ineh Black and Decker 
electric drill connected to it by a leather belt. Rotation 
frequency was controlled by varying the voltage to the drill. 
In order to improve the frequency control at slow speeds 
(below SO RPM) a second polishing-nead spindle was interposed 
as shown in Figure 2. These two configurations then gave a 
speed range for the oscillator unit of from 4 te 1000 RPM. 

A 1/4-ineh plywood plate of &.l4-inches diameter with a 
30-degree cam hump was fabricated. To prevent wear a metal 
strip was cemented to the cam. A similar 10-degree cam was 
alse constructed. The cam, with a 1/2-ineh centering hole, 
was held securely by a nut against the plywood plate on the 
spindle shaft. A compass card oriented to the exact 
escillator degree was secured to this plyweod plate. The 





cam was set at any desired oscillator position by loosening 
the cam-holding nut. 

As shown in Figure 2 a microswitch was positioned adja- 
cent to the cam. The awitch was wired in the normaslily-open 
position to the preset-time relay of a Rediation Instrument 
Development Laboratories, model 206, scaler. The cam was 
positioned by use of the compass card to make contact with 
the awitch at a specified oscillator degree. When the sine 
pattern of the oscillator unit reached this degree the cam 
made contact and the switch was closed. This clesed the 
preset-time relay allowing the scaler to begin to count. 
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When the microswitch broke contact the relay of the scaler 
opened, ceasing the count. 

A Ra@iation Counter Laboratories boron trifloride pro- 
portional counter, model 10503, was insertec midway into the 
assembly through holes cut into the east cover. It in turn 
was connected to the Radiation Imstrument Development 
Laboratories sceler. 

The oscillator unit end selective-count device was in- 
serted from the south end of the assembly until the oscillator 
unit wae at the center. As shown in Figure 4 this placed 
only the oscillator unit itself beneath the assembly; the 
remaining equipment was outside to facilitate chenging the 
can setting. 
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Vi. EXPERIMENTAL PROCEDURE 
Ae Description of Typical Run 


The fuel elements were loaded into the assembly es shown 
in Figure 1. Measurements were taken only in the center of 
the assembly. Holes drillec on the e@ast-cover centerline 
matched assembly holes between the graphite blocks. The 
holes were one foot, two feet, three feet, four feet and 
four feet eleven inches above the base. Nereafter these 
holes will be referred to as holes one through five 
respectively. 

The neutron Getector was inserted midway into the assembly 
and the oscillator unit set at 0 degrees. With the 30- 
degree cam set at a meen position of O deprees, initial 
contact with the microswitch was made at 345 degrees. 

Oseiliator speed was adjusted to the desired RPM by re- 
ferring to a stop watehn and counting revolutions. When de- 
Sired speed was obtained the watch wes etarted and the 
preset-time switch om the scaler was actuated. Runs varied 
in guration from one to five minutes depending on position 
ena oscillator epeed. A run was completed when a designated 
number of counting perioda hed elapsed rather than ending at 
en exact time interval. For example, a run for 3 minutes 
at 16 RPM requires 54 counting periods; when 54 periods had 
@lepsed the preset-count switch was turmmed off even though 





> aimee no .,, ae, : 

















eS ——- Some bool ot 
— «o «Aa mt se mee pee 
p ae ae eke tt 7s 
a ne liad 9 ll emai Siig = a! Soe ae 
a ae ee ee a ee? ee eee le 
a CO LO mm ££: —,, win ae 
-— -—, - ——<— 7c ere «= (ee kee 
7 

Gin imtoo a ma a: aa 
— ee i Ee one om 
CE EE OM nm a 
= ee edie oe! ote lee 

et ml 
em bene eh a me ere 8 et el 
oo. ae — =e or oe me 
at a a a Oe ae me Of 
0 ee Ons GD a + ee ae 
oe 

—) “he eae are we ome 
pom ay hye apt 
ee dss et hee 




















30 


the time may have been a few seconds on either side of 3 
minutes. It was difficult to control the speed exactly and 
this procedure eliminated varying nusbers cf counting periods 
in a single series of mums. 

The counts were recorded from the sealer; the drill 
turned off; the cam reset to a@ mean position of 30 degrees 
and the entire counting procedure repeated. The counting 
procedure continued in 30-degree increments until the entire 
360 degrees had been surveyed yielding one single series of 
runs. 

Speeds (frequency of oscillation) used in the investiga- 
tion were selected at regular intervala on 4 logarithmic 
acale. Runs were attempted throughout the entire rence of 
she equipment. The lowest constant speed attainable was 
i RPM. An upper limit on speed was imposed by the cam and 
microswitch arrangement, for in speeds beyond 160 RPM the 
switch hammered the cam as the switch spring forced the idlcr 
back to the surface of the cam. As rune were made in the 
higher hele positions at high speeds the data beeame unusable. 
This was caused by the exponential attenuation of neutron 
flux up the assembly and a consequently amaller difference 
from the highest to the lowest count rate in a series of 
runs. ‘when this occurred the counting rates when pictted 
against degrees did not yield a satisfactory distribution. 

Criginal data for the various runs completed is shown in 
Tatvies 5 through le in the Appendix. 
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B. Determination of Phase Angle 


The count rates for a single series of runs were plotted 
on graph paper. Standard deviation for each point was also 
plotted to facilitate drawing a representative curve through 
the points. All the plots that showed a trend in distribu- 
tion gave a sinusoidal reault. Figure 5 shows a plot of the 
data for all holes at 2.51 radiens/second, and as mentioned 
previously some unusable data was taken and the basis for 
Geclaring it unugable was large scatter in the plotted 
curves. 

From the general trend of the plot and especially from 
the position of tne hishest and lowest counting rates the 
phase enmsle for that series of runs was determined. To ensure 
that an accurate phase angle was beins obtained 4 10-degree 
¢an was utilised. Five-minute counts were taken every ten 
Gegrees in the region of the maximum counting rate, however 
use of the lO-Gegree cam was restricted to the lower speeds. 

It was sesumed in determining the phase angle that the 
esoillator unit when in the 6-degree position during rota- 
tion and the counting equipment would yleld the highest 
counting rate. Teo validate this assumption an entire series 
of calibration runs with the neutron detector held 1.5 
inches from the oscillater unit in position beneath the 
aasembly was made. Table 13 in the Appendix and Figure 5 
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indicate the results which show that there was no "built in" 


phase angle. 


C. Determination of gf by Amplitude Hatios 


The inverse relaxation length for the thermal-neutron 
flux in the assembly, xy , was determined by amplitude ratios 
for various frequencies of oscillation. The value actually 
obtained was the inverse attenuation length, but as previously 
noted, under certain conditions these values are approximately 
equal. The experimental conditions were such that the above 
statement 1s true, and will be demonstrated in the Discussion 
of Results. 

Uncorrected-ampliitude ratios were used in place of foil 
activities in the equation A = Ce~°*. with the method of 
least squares an initial uncorrected value for y was found 
(2, p. 166). Both end-correction and harmonic-correetion 
factors” were determined and applied to the amplitude ratios 
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with this initial value of y». The origin of the Z-axis was 
taken at the base of the assembly; therefore hole one was at 
@ position of 2 = l2 inches. A second value for 7Y was 
determined with the corrected amplitude ratios. This process 
was repeated until the value of XY used in the correction 
factors agreed with the xy found by using the corrected 
amplitude ratios. Computations showed thet harmonic cor- 
rections were significant only for the three lower holes and 
end corrections for the two upper. | 
It was believed that harmonic corrections were unnecessary; 
therefore a value for XY was determined in a similar manner 


using the end corrections only. 


D. Determination of of by Medilan-flux Levels 


inverse relaxation length was also determined by using 
the median-flux level. Uncorrected median-flux levels were 
used in place of foil activities in the equation A = Ce” v2. 
A procedure was followed as explained in the preceding 
section to obtain @ value for Y . Values were determined 


by using only the end-correction factors. 


EB. Datermination of Amplitude Ratlo 


Basleally, the need of the amplitude ratio was to deter- 
Mine if the amplitude of the neutron sine wave varied with 
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the frequency of the source oscillation. To eliminate any 
“puiltein" variation due to the equipment, the calibration 
runs next to the oscillator unit were used as the basis for 
determining the anplitude ratios. The amplitude ratio was 
defined as the amplitude of the wave obtained for a particu- 
lar position and speed divided by the calibration-run 
amplitude at the same speed. Hoth values of emplitude were 
obtained from the plots of the data, and to standardize the 
results, all amplitudes were placed on a one-minute basis. 

After the amplitude retios were obtained from the data, 
an end-correction factor was applied. From the iterative 
precess explained previously, 2a value for the end-correction 
factor was obtained. No end-correction factors were applied 
to the calibration runs. The end-correction factors were 
Givided into the amplitudes obtained in the hole positions, 
and it was found that factors were significant oniy in the 
top two holes. 

There is a Gefinite difference between the theoretical 
amplitude ratio obtained from Equation 31 and the amplitude 
ratio obtained from the experiment. Theoretical asplitude 
ratio was cGefined as the ratio of amplitude of the observed 
run divided by that of an infinitely slow run, whereas the 
experimentally-determined amplitude ratio is the smplitude 
of the observed run divided by a calibration run. Therefore, 


numerical values of the two ratios wili not agree, but they 
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will agree in showing any change in auplitude ratio with fre- 
quency. 
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¥. Determination of Median-flux Level 


The median-flux level is the average flux at a particu- 
lar hole position and oscillator frequency. From the plots 
of original data, the amplitude and the lowest value of 
count rate found in the graph were obtained. The median-flux 
ievel was then taken as half the amplitude added to the 
lowest count rate. End-correction factors were applied to 


the amplitudes for the top two holes. 
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Vil. RESULTS 


The original data for all runs completed including the 
calibration is given in Tables 5 through 13 in the Appendix. 
¥isure 5S shows a representative plot of original data at a 
frequency of 2.51 radians/second for all hole positions and 
calibration. 

Phase angle was determined for each run and the results 
are as given in Table 1. All the calibration runs proved 
thet there is no “built in” phase angle in the equipment. No 
phase angle is observed in any hole position above 0,63 
radians/second. Approximately ea 20-degree phase angle is 
found to exist in ali five hole positions at 0.63 radians/ 
second. For the only hele position investigated, number two, 
a 30-degree phase angle is found at 0.42 radians/second. 
Verification of the phase angies in hole two up through 4.19 
radians/second was obtained by use of the 10~-degree cam. No 
change of phase angle is observed in proceeding to higher 
positions while taking measurements at the same oscillator 
frequency .« 

Amplitude ratios were obtained and the results are as 
given in Table 2. Figure 6 is a plot of measured amplitude 
ratio versus csclllator frequency for hole positions one, 
three and five. Amplitude ratios do not show a definite 
trend with frequency, but do remain approximately constant 
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Table 1. Phase angle (degrees) 





Hole position 


Frequency ee ae etn a renee nemesis mer amen SY ATI ae 
radians/second One Two Two" Three Four Five tion 
0.42 — » -~_ «= © 
6.63 2Q 20 ao 20 20 20 0 
1.05 Q GO G 8) o 0 QO 
1.65 G G 0 0 Q Q QO 
2.51 Q 0 $) G Q ) 0 
4.19 Q Q Oo 6 9) 0 QO 
6.28 oO G oa 0 — — ©: 
7-33 0 GC «= oo -~ = *) 
9.42 $$ + == ._ = o- Oo 
12.55 5 ¢ — — dx 0 
18.55 i- - —_—-— = O 





“guns accomplished with 1lO-dgree cam. 


Table 2. Amplitude ratio 





Hole position 


Frequency 
radians /second One Two Three Four Five 
0.63 0.0710 0.0330 0.0125 0.00610 6.00418 
1.05 0.0076 0.0440 0.0139 9.00663 0.00279 
1.68 0.0786 0.02438 0.0083 0.00450 0.006216 
2.51 G.0510 6.0150 6.0090 0.00492 0.00274 
4.19 0.0696 040315 6.0109 6.00454 0.00252 
6,28 0.0550 0.0255 0.0083 — ~—) 
7033 0.0536 — or = cow on on 
93 2 0.0673 0.0154 — — — 
12.55 0.0673 0.0218 -- a =o 


18.55 0.9715 — saadaeal =o wan 
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Figure 6. Amplitude ratio vs. frequency 
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at each hole position. A high coefficient of variation is 
evident for values above hole one. 

Amplitude ratios were also plotted against vertical 
distance in the essembly. The five frequencies investigated 
in all five hole positions are represented on this plot, 
Figure 7. The line slopes drawn correspond to the x values 
obteined by the method of least squares for the set of 
points used. 

& median-flux level for ¢ach hole position and frequency 
was determined. Table 3 gives the results of this computa~ 
tion, which also includes the calibration runs. To demon- 
strate the change in median-flux level with oscillator fre- 
quency 2 plot of these two varlabies was constructed for the 
first three hole positions end calibration run. Figure & 
Shows that the slopes on logarithmic paper of tnese lines 
for the four positions are equal. 

Table 4 gives the various values of Y obtained by using 
the amplitude ratios. Uncorrected values gave an sverage 
ef 0.0710 inchee™’; end-and-harmonic-corrected values gave 
an average of 0.0631 inches™*, and end-corrected values gave 
an average of 0.0702 inches”*. End-corrected y values were 
plotted versus frequency in Figure 9. Values for ry are 
approximately constant with frequency. Large coefficients 
of variation are noted for all the Y values plotted in 
Figure 9. 
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Table 3. Median-flux level (total counts) 





Frequenc eee Calibra- 
radians/” _~ Two Thre + ~, = 


0.63 11,534 5,410 3 » 935 2,105 1,150 13 585 


1.05 11,275 5,635 39 (30 2,115 1,265 135425 
1.03 11,350 B22 (0 »473 1,996 2085 3,245 
2.5, 11,150 5,315 3,590 2,098 1,183 12,947 
4.19 10,857 5,11 33450 2,160 1,098 12,739 
6.23 3 &8O 5,200 3500 ae —— 12,536 
7.33 10,634 ~ == ~~ -- ~- 12,495 
9:42 10,600 5,078 - —~ —- 185339 
12:55 10,490 4/830 - - ~ 
18.85 10,400 —_ — eal —_ 11,935 


Table 4. Inverse-relaxation length (inches )7+ 





Method 

of “4 
panrabe- Frequency-radians/second 

tion 0.63 1.05 1.68 2.52 4.19 Average 
No cor- 
rections 0.0642 6.0714 0.0764 0.0690 0.0739 0.0710 
End cor- 
rections 0.0640 6.0704 06.0759 0.0676 %&.0729 4.0702 
End and 
harmonic , 
cor- 0.0582 6.0636 ¢.0690 6.0626 6.0664 0.0631 
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Figure 8. Median-flux level vs. frequency 
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A value for xy of 0.C0G5l1 ineches™* was obtained by using 
the method of median-flux levels. 
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VIII. DISCUSSION OF RESULTS 


It was the original belief that a phase angle would be 
observed at any oscillator frequency above 0.6 radilans/second 
(13, p» 33). After a series of runs had been completed, this 
was found not to be true. Previously, the values of k and 
head not been kmown for the assembly used, and had to be 
esoumed. This led to the erroneous prediction. 

the theoretical enslysis was performed using one everage 
ail show the 





group of delayed neutrons. Figures 10 al 
theoretical phase angle and amplitude ratio respectively eas 
Getermined from Equation 31. <A value for k of 0.53 was used “ 
With f equal to 6.87 (9, p. 56), and B* equal to 76 x 

107° om*® (9, p. 57d) a value of £ was determined to be 
6.00098 seconds from the equations 


et. “fis Lim 
Li L“(1*f) and 1 + Lp" 


and an assumed value of /, = 0.001 seconds. The total 
fraction G of the delayed neutrons was set equal to 0.0064 
(7) with a corresponding value of . equal to 0.08 seconds™*. 
It can readily be seen from Pigure 10 that a phase angie 
would not be observed in the assembly until en cecillator 
frequency of 40 radians/second is reached. A phese angie 





“Beck, D. M., Ames, Iowa. Data from partially completed 
M.S. Thesis. Private conmmunication. 1953. 
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does not appear for this investigation until very high fre- 
quencies because the value cf k is low. If kK were 0.9 a 
phase angle would be observed at a frequency of 20 radians/ 
second. The overail shape of the theoretical curve remains 
the same with varying values of k. 

Runs were completed to the lower and upper frequency 
iimits of the oscillator and counting equipment. The upper 
limit (19 radians/second) is short of the required frequency 
for an observable phase angle. In complete agreement with 
theory, no phase angle was observed in any of the hole posi- 
tions above 0.63 radians/second. Use of the 10-degree cam 
verified these results. 

Readings were taken in all hole positions for a frequency 
of 0.63 radians/second. A consistent phase angle of approxi- 
mately 20 degrees is found. At 6.42 radians/second the phase 
engle in hole two is approximately 30 degrees. The theory 
advanced does not predict a phase angie at these low valuee 
of frequency. Further investigations were not conducted into 
this low frequency range since the equipment did not provide 
adequate speed control below 4 RPM, 

Theory makes no distinction of position in relation to 
phase angle; therefcre the phase angle should be identical 
in ali areas of the assembly. No change in phase angle with 
hole position was observed during the investigation. Surveys 
were taken only up the center ef the assembly se the 


i Om — le hie! ele a) ere eee cin: 
he Me S Neer ew + eee Bee 
= © } eet © te frre of Slee Ghee cue, 
ele ae cee ee |e 
a’ = cleus cus eae oe 

—7 a A ed et «ct lois Gate comet 
—=— Geek Bm) ak tlle ot DS ota 
ey! ee el @ ee. emer en Wt A 
Te ee th Fe ee ee 
dew Oh ee tS ak betes coe ee ee ee 
ed ee! ee cen Ola 1 ml 
fe) mall thE omg 

OO) 6 6) ee sien Met ho ff cau #78) aot 

ee ee 
ee eee ee ee 
me ek ye 6 oe ola we 
lee mel me em ee Ou Om 

en iw Tete lem ot Ae ld eg 
Ory we A Ame a ee TE! mk 
et ome eee boom 
RR dE Leg FS rT Oe pe 
a a Cte eet omens me) ee jolr « 
Ae ee ee Ce a 
oo en OF Wy Lene eet ge Lm meee 

















a 
S 
, 








qe 


Ra 












52 


proof of the constancy of phase angle ia Limited to this con- 
dition. With the method of source forcing used it is believed 
that some large variations would be found at the edge cof the 
assembly near the base. Simply noting that the pattern 
position facing these extreme positions will vary as much 

as 4&5 degrees from the recorded position indicates the 
prysical reason for the nonconfermance. 

Ag shown in Figure 5, the curve of the neutron flux 
measured at the various hole positions is shown to be 
Sinusoidal. Scatter of pointa increases as higher hole 
positicns are investigated. The range of frequency for 
usable data narrows as measurements are taken higher in the 
assembly. Statistical variations, and inaccuracies in the 
Getector and scaler distort the counts received. Due to the 
already eamall amplitude of the sine wave, any distortions 
cause a scattering of points that make the readings unusable. 
Thies situation is agsravated by higher frequencies. As the 
osclliator frequency increases, the neutron waves are propa- 
gated quite close together. Any disturbance or 
irregularity in the entire assembly or counting system will 
cause an incorrect measurement to be taken. For example, a 
measurement is to be taken at a specified oscillator-unit 
Gegree. A small disturbance or irregularity will cause the 
néutron wave to be delayed in reaching the detector. At a 
high frequency this delay is significant since the portion 
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of the sine wave actually surveyed may be 30 to 40 degrees 
removed from the portion desired. 

As shown in Figcure 5, the standard deviation at hole one 
4s reasonable and acceptable. However as readings are taken 
at higher positions in the assembly the coefficient of varia~ 
tion caused by standard deviation increases even thoush 
longer counting times are allowed. This is caused by lower 
count rates due to the exponential attenuation of the thermal~- 
neutron flux. An average coefficient of variation for all 
the runs completed renges from cne to three percent in going 
from hole one to five. It is seen then that to obtain 
eoefficients of variation of less than two percent at the 
higher hole positions, very long counting times are needed, 
in the order of 10 minutes at hole five. 

Amplitude ratios were determined and are given in Table 
2. These amplitude ratica for holes one, three and five were 
plotted against frequency in Pigure 6. A least squares analysis 
was applied to the data for holes three and five, and it was 
found thet the line joining hole-three data had a slope of 
“0.104 40.225, and the line joining hole-five date had a slope 
of -0.231 40.344. This indieates that a horizontal line 
representing a constant value of amplitude ratio with frequency 
lies easily within the standard deviation for the slopes found 
fer both holes. It 4a believed that the amplitude-ratio values 
are approximately constant as predicted by Figure 11 and that 
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the slopes shown on Figure 6 are caused by the high standard 
Geviation rather than varlance of emplitude ratio with fre- 
quency « 

Figure 11 has the overall shape that the amplitude ratio 
versus frequency would have for any value of k. For ak 
greater than 0.53, an amplitude ratio less than one would be 
noted at frequencies less than indicated on Fisure 11. 

Seatter of the pointa in Figure 6 is due to a very high 
coefficient of variation. As was mentioned previously, some 
seatter of points ise present at all positions especlally in 
the higher hele positions at higher frequencies. Error occurs 
wnen fairing in the curves to f1t the data. A very signifi-e 
cant error cecurs when the standard deviation of the points 
is large compared to the amplitude. For example, at hole 
three a standard deviation of approximately 66 coun 
found; however an amplitude of only 300 counts exists. 
Coefficient of variation caused by standard deviation of the 
amplitude is then approximately £O percent. The condition is 
exaggerated at hisher hole positions such that a 36 to 42 





percent coefficient of verlation in hole five is noted. 
Significant errors appear for the amplitude ratios above 
hole one. Considerably lonser counting times are needed at 
all higher hole positions to give satisfactory results. 

A counting time of approximetely 16 minutes would give a 13 


percent coefficient of variation for hole four. 
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Even though the coefficlents of variation for the ampli- 
tude ratios are quite high, the results are believed to be 
sufficiently accurate and repetitive to support the conclu- 
sions drawn. 

Values of the median-flux level are given in Table 3 and 
the plot of level versus frequency is shown as Figure &. 
Deviation was so slight in holes four and five that they were 
net included in the figure. 

The transfer-function development does not predict a 
decrease in the median flux, yet all four curves show that 
mecian-flux level does Gecline as frequency increases. Slopes 
on the logarithmic plot (Figure &) of all four lines are 
found to be equal. If each line ls represented by the equa- 
tion = Ki)", where K is a constant end n is the slope of 
the line, then the rate of neutron-flux deciine with fre- 


quency is 

ag Kn) ® ng 

dW Ge) ce , 
The velative rate of deciine is 


ad n 
~~. a. 


enmenncaaitionnpan * CD e 





With identical values of n, and the same values of (, each 
hole position has the game relative rate of decline. The 
Gecline is therefore not a function ef the assembly. 
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The decline is attributed to either the method of source 
forcing or to the counting equipment. No logical reason can 
be advanced for the unit itself to cause the decline, but 
the counting equipment could easily have causecé this un- 
predicted result. The number of activations of the preset- 
count relay is in proportion to the frequency. A small relay 
lag, would then become important at higher frequencies and 
Significantly reduce the amount of time the scaler was count- 
ing. Furthermore, investigation shows that the reduction of 
thermal neutrons and fast neutrons counted is in proportion 
to their total counts. The deviation from theory is there- 
fore said to be a mechsnical irregularity of the counting 
equipment rather than a function of the assembly. 

Amplitude ratios were not affected by the decline of 
counts with frequency since the calibration runs were used 
as the denominator of the ratio, 

When plotted on semi-log paper, amplitude ratio versus 
vertical distance gives a reasonably straight line. Figure 7 
shows the results for frequencies from 0.63 to 4.19 radians/ 
second. As noted previously, the large coefficlents of 
variation of data for the higher hole positions result in 
wome scatter of the amplitude ratios. It ls seen, however, 
that the points nearly all lie close to the straight lines 
indicating that the points can be considered reliable. 
Standard deviation of the amplitude ratics is indicated on 
the figure. 
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Values for inverse relaxation length, XY , were computed 
for the frequencies. These values obtained by an iterative 
process are listed in Table 4, and it is found that values 
range from 0.0582 to 0.0759 inches™*. A value for Y¥ was 
fourid without using elther the end or harmonic corrections; 
the average value was ©.0710 inches7*. When both corrections 
were applied, an average value of 0.0631 inches~* was found. 
It is beiieved however that the higher harmonics have been 
greatly iessened with this particular arrangement of sources 
in the oscillator. The sources are placed in a line and are 
surrounded by three inches of paraffin which thermalizes 
many of the neutrons produced. The paraffin, lined sources, 
and two water-filled tanks in the outermost pedestal spaces 
help to scatter the neutrons more uniformly beneath the base, 
so that a plane source of thermal neutrons is approximated. 
Dopchie, Leonard, Neve de Mevergnies and Tavernier (3) 
showed in their work with an exponential assembly that dis- 
tributing four sources in a square form reduced the harmonic- 
correction term to 5 percent or one-fourth the correction 
needed when the sources were placed together. 

With application of just the end corrections an average 
value for ¥ of 0.0702 inches~* was determined. Hayes (8, 

p. 27) found by the method of foil irradiation that the 
assembly in this configuration had a Y of 0.0705 inches™:. 


Harmonic-correction factors would appear then not to be 
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needed, or if used the full correction should not be applied. 
Line slopes on Figure 7 are the XY values determined by using 
only the end corrections. 

Values of inverse relaxation length obtained by using 
end corrections only were plotted with frequency on Figure 9. 
A standard deviation for yY was obtained by using the method 
of propagation of precision indexes as presented by Werthing 
and Geffner (17, p. 209) and are indicated on the figure. 
Large coefficients of variation are again noted for 2. 
This is due to the high ccefficilentsa of variation of the 
anplitude raties, especially at the higher hole positions. 
By using longer counting times the relative standard devia- 
tion of the amplitude ratios would be greatly reduced and 
statistically better values of og would follow. For example, 
if an average relative standard deviation of 5 percent were 
found for the amplitude ratios, then the coefficient of 
variation for y in this experiment would decrease by 2 
factor of approximately six. 

it is seen from Figure 3 that there is no trend in x 
with increasing frequency. A constant value of XY should be 
Found, but due to the high standard deviations values differ- 
ing slightly from the average of 0.0702 inches™* are noted. 

To substantiate the values of XY determined by using 
amplitude ratiog an average value of Y was found by using 
the median-flux levels. As shown on Figure l2 the values of 


Ve 


xii ieee eet they av Wee Uw eee 
ram , ke he T cue) et ode ted 
ae ier ate a) 
oe amie 
we fetes © eegues Mmuselse, owe 
—— erp) LN ends rn ot ater! ep 
swu= oc sabe of benlevm oa | Wey eabvelre Cee, 
oe ence 
rol oer ae emer ee See (CS wy (WT! el 
~~ oe oe ae ara oe ee 
eab @ becinies US couiint Vee Ge ae ot Be Ot 
ee ee 
ann mmol vo eee pie B 
oY TF he eee suiiew «ty Is ad 
Aas A OU Ge) me ee ee ee ee 
ew eres, = a Ue ee woh) epee a BS 
& penicltws oo aame eddie: stetiians om wt comet 
. owe ee eee ee we) eee 
ea At ys mare 
ee oF ad ee ee et ee oe a 
Call wire me + Ca Gelewewrels MOET 
a pon neni a Rey BT 
— as 
et See  Qeeeye aa arr, 
; ea er ee ce 
—— & eoeeve * | 
oc ee te) lee eee a la oe tIC® 
a 
mr we En ne ote 









MEDIAN-~FLUX LEVEL ~- COUNTS/MINUTE 


59 


6000 r 
al 0 = 0.63 RADIANS/SECOND 
e = 1.05 RADIANS/SECOND 
4 = 1.68 RADIANS/SECOND 
4000 a = 2.51 RADIANS/SECOND 
Vv = 4.19 RADIANS/SECOND 
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Figure le. Median-flux level vs. vertical distance 
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the median flux are grouped close together. An average value 
was obtained from each group and a value of oS computed fron 
them by using end corrections only. Standard deviation of 
the average points is approximately the sige of the symbols, 
therefore the value cf XY obtained is reliable. From Figure 
12 s Value of ¥Y was found to be 0,068] inches™~* in compari- 
son to the average value of 0,0702 inches”* computed by using 
amplitude ratios. Agreement within 3 percent between the 

two methoda indicates the amplitude-ratio method to be 
acceptabie. 

The validity of the statements concerming the attenua- 
tion length end 1/y wes investigated. Values used in the 
computation were T = 350 om® (16, p. 331), Ky = k(1 - 6) = 
0,526, L® = 325 em*, £, = 0.001 seconds and GW = 20 
radians/second. A value was determined from Equation 32 for 
attenuation length and found to be 14.3 inehes. The inverse 
of this number is 9.0676 inches~*+ which agrees favorably 
with the values of C.0702 inches~* and 0.06681 inehes~* ob- 
tained by the methods of amplitude ratioa and median-flux 
Levels respectively. 

Frequencies utilized in the experiment were small enough 
to cause the attenuation length to be approximately equal 
to the relaxation length. From Equation 32 it is noted that 
a frequency of 20 radians/second produces negligibe change 
in the value of attenuation length as compared to very low 
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frequencies. Hence from the previous computations and dis- 
cussion it appears that the attenuation length and conse-~ 
quently the complex buckling for the frequencies investigated 
are approximately the same as relaxation length and material 
buckling. Therefore 1% would be proper in this case to 
utilize the values of attenuation length for determining Y . 
Even with the combination of varying ana steady neutron flux, 
eomputetion of XY by amplitude ratio and median-flux level 

is feasible since the frequency never is large enough to 
produce significant changes in YY . 

For assemblies with much larger values of Kos frequericies 
used in this experiment would produce changes in the attenus- 
tion length and give varying values for the complex buckling, 
therefore the inverse attenuation length would not be equal 
to the inverse relaxation length. 
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IX. CONCLUSIONS 


A theoretical and experimental investigation into neutron 
aiffusion in the subcritical assembly with sinusoidal source 
oscillation was conducted at Iowa State College. An oscillia- 
tor unit produced the sinusoidal source osciilation of the 
thermal neutrons, while a selective-count device allowed 
discrete sections of the sine wave to te cbserved and counted. 
Neutron-diffusion information was extracted from plots of 
these counts. Conclusions drawn from the information are that: 

1. No phase angle is observed for oscillator frequencies 
between 0.63 and 15.85 radians/second. This is in agreement 
with theory. 

2. An unpredicted phase angle of 20 degrees is observed 
at a frequency of 0.42 radians/second and an angle of 30 
degrees at 0.63 radians/second. 

3. There is no phase~angle change at higher hole posi- 
tions on the centerline of the assembly at the same oscillator 
frequency, which ranged from 0.63 radians/second to 6.28 
radians/second. This is in agreement with theory. 

4, Amplitude ratios for a given hole are approximately 
constant for varying frequencies from 0.63 radians/second 
to 6.25 radians/second. This is in agreement with theory. 

5. The median-flux level decreases as frequency in- 
creases. Mechanical operation of the counting equipment 
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causes this change rather than any assembly peculiarities or 
characteristics. A change in ievel is not predicted by 
theory. 

6. The osciliator unit thermalizes enough neutrons to 
give satisfactory readings up through hole five up to 4.19 
radians/second. Any frequency readings above 12.55 radians/ 
second are confined to nole one. 

7. An inverse relaxation length of 0.0702 inchea~* 
was Getearmined by using amplitude ratios. 

&. An inverse relaxation length of 0.0651 inches7™* was 
determined by using median-flux levels. 

9. The atterwation length of the neutron wave is equal 
to a constant value of l/y for the assembly and frequencies 
investigated. 

io. Longer count times are needed to reduce coefficients 
ef variation in the ampiitude ratios. in ail other respects 
the experimental equipment and preeedure will yield satis~- 
factory reaults within certain frequency limits for phase 
angle, amplitude ratio and inverse relaxation length. 
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X. SUGGESTIONS FOR FURTHER STUDY 


From a review of the literature it is obvious that neutron 
diffusion in a subcritical assembly with sinuscidal source 
oscillation has not been fully investigated. 

yhe equipment used in this experiment can be improved in 
several ways. Another epeed reducer could be introduced to 
S&llow speeds below 4 RPM to be utilized. The can could be 
replaced by a phototube device that would trigger the scaler 
at the higher frequencies. An automatic revolution counter 
would be necessary to measure the frequency at the higher 
speeds. A gating circult would be more desirable to inter- 
rupt pulses going to the sealer than using the preset-time 
relay which 1s burdened at high frequencies. Still another 
possibility might be to use an ionisation chamber connected 
to a sensitive electrometer. The wave produced by a re- 
gorder connected to the electrometer could then te analyzed 
for data desired. 

In this experiment, only the center of the assembly was 
investigated, Measurements could be taken throughout the 
eseembly to see if the thermal-neutron Slux behaved in the 
by a steady source. Further 
theoretical and experimental work with neutron weves could be 
accomplished with an investigation into the change in prompt 
buckling by souree forcing. Speeds considerably below and 
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above those investigated here cculd aiso be utilized to ex- 
pand the experimental verification of the theoretical transfer 
function. The above suggestions could be expanded greatly 

by changing the lattice in the assembly and repeating the 
particular inveatigation. Voids and absorbers could also 

be placed in the assembly to observe the effect of them upon 


the neutron wave, 
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Table 5. Original data at 0.63 radiana/second (total 
numbe ts ) 








r of coun 
Oscillater Hole position 
position 
One THO Three Four Five 
degress (3 min.) (3 min.) (3 min.) (4 min.) (5 min.) 
0 17 2900 8,623 3,794 2,254 1,234 
2° 15,531 &,71il 3,705 2 9220 1,257 
15,310 § ,600 3,660 @ 9170 1,202 
Go 18,056 G2e4 3,040 2,105 +229? 
150 16,306 7,Ck7 33522 2,031 1,074 
1360 16,026 72708 3,374 2921) 1,042 
210 15,99 7752 3368 1, 1,100 
240 16,24 757&5 3,484 2,026 1,077 
270 17,093 7999S 3 491 2,036 1 ,OS6 
330 17,747 & O56 3,685 2,118 2204 





fable 6. CGriginel data at 1.05 radians/second (total 
rmunser of counts) 


a 





CRE PRL 





Oselliator Hole poaition 
— one Two Three Four Five 
& rs 

degrees (2 min.) (2 min.) (3 min.) (4 min.) (5 min.) 
e) 11,950 6,342 3,96 2,248 1,297 
30 12,927 6,036 3,937 2,293 1,370 
60 12,681 5,f2l 3716 2,165 1,309 
90 11,472 54386 3572 2,215 1,282 
120 42,037 Se s5L 3319 2,193 1,201 
pf 10,749 59230 at 2,019 L209 
i rag * 5 ye2e3 3,498 1,968 1,200 
210 10, 5 ores 30982 1,934 1,242 
240 11,207 32233 3506 2,034 1,236 
276 11,777 559207 $992) 2,331 1,255 
306 11,210 e478 3% 2; 1,229 
330 11,373 § TBS 3,579 @ 3239 2,249 
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Table 7. Original data 2 radians/second (total 


Ti 


number of counts 


OseLllator Hole position 
position 
One Two Three Four Five 
degrees (2 min.) (2 min.) (3 min.) (4 min.) (5 min.) 
a Ce 
: 2 P a tf 2 
90 11,472 5,25 3+pe 2,067 1,099 
120 12,065 5,282 3,444 1,959 1,11 
1 10,59 peli’ 3,460 1,932 1,08 
1 10,625 982 3,310 1,92 1,076 
219 10,572 5» 000 3%39 1,694 1,035 
240 10,946 5 214 3,369 1,951 1,025 
270 11,291 5,191 39925 1,921 1,03 
360 11,665 59189 39479 128 1,416 
330 12,257 55234 3,499 1,902 1,131 


number of counts 


Table &. Original data — radians/second (total 





Osciliator 
position 


degrees 
30 
60 
GO 
is 
160 
210 


240 
2fo 


300 
339 





Hole position 


2,196 


2,177 
2,045 
2,034 
2 O00 
2,002 
1,975 
2,092 
25150 
2,136 


One Two Three Four Five 
(2 min.) (2 min.) (3 min.) (4 min.) (5 min.) 


1,235 
2 





eg lm A 
— “se 


| 





Table 9. Ordginal data So radians/second (total 


72 


number of counts 





Oscillator 
position 


degrees 


0 
eo 
90 


129 


Kole position 


Table 10. Original data at 6.28 redians/second (total 
number of counts ) 





Oscillator Hole positicn 
position 
One Two Three 
Gegrees (l min.) (1 min.) (3 min.) 
50 Bros «eka «08 
29 {US ac 23D 
GO O80 3,256 & O12 
90 aes 3,120 3,908 
120 5,641 3,013 3,857 
150 5933702959 3,823 
186 5,272 2,981 3,520 
219 5,334 3 O50 3,830 
240 2521 3,129 3,067 
27o —— 34130 32897 
300 5s (20 gOS tees 
330 Os 155 3,245 9056 
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Table 11. Original data at hole position one (total number 





of counts ) 
Oscillator ape 
position “— 
re 
deprees 3 _ (2 on (2 min.) a “it 
9) 269347 11,261 5,500 5,357 
30 16,47! oat oaks 5,452 e834 
GO 15,9 9 3 5414 52359 
15,005 102344 59 354% 5,161 
120 15,246 10,458 92007 h 985 
130 15,288 9,962 O94 4d 
210 25,320 10,140 4,929 4,916 
eho 15,932 10,729 BLS 5,241 
270 16,402 10,970 5,218 5 326 
CO 16,390 11,133 $5508 § 516 
330 9955 11,190 % 5605 5308 


Table 12. Original data at hole position two (total number 





of counts) 
Oscillator 
position Radians/second 

7.33 9.42" 9.42 12.55 

degrees (3 min.)  —- (2 mans) 2 min.) (1 wld. 
0 7, 5a4 eS 137 5,152 & 951 
30 7,089 5058 55197 4906 
90 7382 prone 5,056 be97S 
420 7,290 9303 5116 f 687 
150 7,396 4,975 5,111 % 688 
AcG T9577 4 029 4 ¥ 906 4,573 
216 75388 4,870 ‘i ,o72 & 751 

246 73503 4,975 9 hO3 475 
270 7,677 5,147 5,239 GUS 
300 75551 52332 5,201 4 938 
330 73539 5,232 5,270 5,063 

Spiret run. 


“second run. 
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Table 13. Original data for calibration (counts/ninute ) 
Cselllator 
position Redians/second 
degrees 0,42 0,63 1,05 1.68 251 #19 
9) 16,658 193,215 16,956 18,700 18,216 18,003 
30 17,374 16,291 16,335 18,757 16,269 17,513 
60 10,174 17,270 17,366 17,475 17,377 17,016 
go 14,729 16,850 p bh 15,565 15,615 14,697 
129 13,227 13,609 13,153 12,324 12,576 12,442 
150 95155 9,543 9,701 Ps 9 9453 6,957 
150 5,049 7,947 7, [o(90 Js677 75800 
210 8,573 &,378 68,651 &,207 372 «©6169 
240 10,693 11,497 11,333 10,909 10,955 10,726 
270 wo 14,432 14,657 14,104 13,881 13,592 
330 ties 1 9352 18,198 L7,777 17,763 17,149 
degrees 6.23 7.33 9.42 12.55 153.85 
O 17,895 17,897 175353 17,542 17,162 
30 16,954 17,761 104952 17,125 eae tt 
g0 13,244 14,208 12,912 233,746 13,523 
1.50 7:783 6,505 7,407 8,419 7,923 
Loo fohTe 7,995 3939 6,955 6,799 
210 S,429 7,753 8,401 &,063 8,144 
240 11,594 10,6059 11,391 10,1? 10,844 
300 16,032 16,303 16,077 15,862 15,3389 
539 17,049 17,642 17,3038 16,967 16,353 





et hee een La ee oe — Ps Voz? 


























